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We compute the vacuum polarisation correction to the binding energy of 

nuclear matter in the Walecka model using a nonperturbative approach. We 
first study such a contribution as arising from a ground state structure with 
baryon-antibaryon condensates. This yields the same results as obtained 
through the relativistic Hartree approximation of summing tadpole diagrams 
for the baryon propagator. Such a vacuum is then generalized to include 
quantum effects from meson fields through scalar-meson condensates. The 
method is applied to study properties of nuclear matter and leads to a softer 
equation of state giving a lower value of the incompressibility than would be 
reached without quantum effects. The density dependent effective sigma mass 
is also calculated including such vacuum polarisation effects. 
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I. INTRODUCTION 

Quantum Hadrodynamics (QHD) is a general framework for the nuclear many-body 
problem |lj-[|. It is a renormalisable relativistic quantum field theory using hadronic degrees 
of freedom and has quite successfully described the properties of nuclear matter and finite 
nuclei. In the Walecka model (QHD-I) with nucleons interacting with scalar (a) and vector 
(u) mesons, it has been shown in the mean-field approximation that the saturation density 
and binding energy of nuclear matter may be fitted by adjusting the scalar and vector 
couplings M. This was first done by neglecting the Dirac sea and is called the no-sea 
approximation. In this approximation, several groups have investigated the effects of scalar 
self-interactions in nuclear matter || and finite nuclei || using a mean-field approach. 

To include the sea effects, one does a self-consistent sum of tadpole diagrams for the 
baryon propagator 0. This defines the relativistic Hartree approximation. There have also 
been calculations including corrections to the binding energy up to two-loops || , which are 
seen to be rather large as compared to the one-loop results. However, it is seen that using 
phenomenological monopole form factors to account for the composite nature of the nucleons, 
such contribution is reduced substantially || so that it is smaller than the one- loop result. 
Recently, form factors have been introduced as a cure to the unphysical modes, the so called 
Landau poles [ IU| , which one encounters while calculating the meson propagator as modified 



by the interacting baryon propagator of relativistic Hartree approximation. There have been 
also attempts to calculate the form factors by vertex corrections [|ITJ. However, without 



inclusion of such form factors the mean-field theory is not stable against a perturbative loop 
expansion. This might be because the couplings involved here are too large (of order of 
10) and the theory is not asymptotically free. Hence nonperturbative techniques need to 
be developed to consider nuclear many-body problems. The present work is a step in that 
direction including vacuum polarisation effects. 

The approximation scheme here uses a squeezed coherent type of construction for the 
ground state [[L^HM] which amounts to an explicit vacuum realignment. The input here is 



equal-time quantum algebra for the field operators with a variational ansatz for the vacuum 



structure and does not use any perturbative expansion or Feynman diagrams. We have 
earlier seen that this correctly yields the results of the Gross-Neveu model [Q as obtained 
by summing an infinite series of one-loop diagrams. We have also seen that it reproduces 
the gap equation in an effective QCD Hamiltonian [[D| as obtained through the solution of 
the Schwinger-Dyson equations for the effective quark propagator. We here apply such a 
nonperturbative method to study the quantum vacuum in nuclear matter. 

We organise the paper as follows. In section 2, we study the vacuum polarisation ef- 
fects in nuclear matter as simulated through a vacuum realignment with baryon-antibaryon 
condensates. The condensate function is determined through a minimisation of the thermo- 
dynamic potential. The properties of nuclear matter as arising from such a vacuum are then 
studied and are seen to become identical to those obtained through the relativistic Hartree 
approximation. In section 3, we generalise such a vacuum state to include sigma conden- 
sates also, which are favoured with a quartic term in the sigma field in the Lagrangian. The 
quartic coupling is chosen to be positive which is necessary to consider vacuum polarisation 
effects from the sigma field. We also calculate the effective sigma mass arising through 
such quantum corrections as a function of density. The coupling here is chosen to give the 
value for the incompressibilty of nuclear matter in the correct range. Finally, in section 4, 
we summarise the results obtained through our nonperturbative approach and present an 
outlook. 

II. VACUUM WITH BARYON AND ANTIBARYON CONDENSATES 

We start with the Lagrangian density for the linear Walecka model given as 

C = ^(iffy -M-g a a- gsfwM 

+ l -d»od,a - l -m 2 a a 2 + ^m>^ - \u^u^ (1) 

with 

£*V = d tl u v - d v u^. (2) 



In the above, if), a, and uo^ are the fields for the nucleon, a-, and cu-mesons with masses M, 
m a , and m u , respectively. 

We use the mean-field approximation for the meson fields and retain the quantum nature 



of the fermion fields MM. This amounts to taking meson fields as constant classical fields 



with translational invariance for nuclear matter. Thus we shall replace 

9°° -> (go°) = go°o (3a) 

9u^n -» (9u^(i) = 9wU^ = QojUJo (3b) 

where (• • •) denotes the expectation value in nuclear matter and we have retained the zeroth 
component for the vector field to have nonzero expectation value. 
The Hamiltonian density can then be written as 

TC = TCn + Ti-a + 7i,u (4) 

with 

H N = 4>\-ia-\7 + (3M)ip + g a aifjif} (5a) 

H a = ^m 2 a a 2 (5b) 

Hu = g<j*wl)H - 2 m !^o ( 5c ) 

The equal-time quantization condition for the nucleons is given as 

[ip a (x, t), ipl(y, t)]+ = 5 af3 5(x - y), (6) 

where a and (3 refer to the spin indices. We may now write down the field expansion for the 
nucleon field ip at time t = as given by |L6] 



^{x) = — ^ J [u r (k)c Ir (k) + V s {-k)c Is (-k)} e lls dk, (7) 

with ci r and ci. s as the baryon annihilation and antibaryon creation operators with spins r 
and s, respectively. In the above, U r and V s are given by 



U r {k) 



( cos^ 



\ 



a ■ k sin ^^ / 



u Ir ; V s (-k) 



( -ff.k s [ n m\ 



V 



cos 



x(k) 



Vis 



(8) 



/ 



For free massive fields cosx(k) = M/e(k) and smx(k) = \k\/e(k), with e(k) = y k 2 + M 2 . 

The above are consistent with the equal-time anticommutator algebra for the operators 
c and c as given by 



c Ir (k),c\ s (k')} + = S rs S(k - k') = [c Ir (k),c\ s (k')] + . 



st (t'\ 



(9) 



The perturbative vacuum, say | vac), is defined through cj r (k) \ vac) = and c\ r (k) \ vac) = 
0. 

To include the vacuum-polarisation effects, we shall now consider a trial state with 
baryon-antibaryon condensates. We thus explicitly take the ansatz for the above state as 



\vac') = exp / dk f(k) c\ r (k) a rs cj s (—k) — h.c. \vac) 
= Up\vac), 



(10) 



di(Jfe) 
di(-k) 



^ I Cj(jfc) ^ 



(11) 



Here a rs = u\ r (a ■ k)vi s and f{k) is a trial function associated with baryon-antibaryon con- 
densates. We note that with the above transformation the operators corresponding to | vac') 
are related to the operators corresponding to \vac) through the Bogoliubov transformation 

cos f(k) — a ■ k sin f(k) 
,a-ksmf(k) cos f(k) J \ci(—k) 

for the nucleon. 

We then use the method of thermofield dynamics [T7J developed by Umezawa to con- 
struct the ground state for nuclear matter. We generalise the state with baryon-antibaryon 
condensates as given by fllQf) to finite temperature and density as Jl3 



/ 



\F(p)) = U(P)\vac f ) = U([3)U F \vac). 



(12) 



The temperature-dependent unitary operator U{(3) is given as JT7 



U(/3)=exp(B^P)-B(P)l 



(13) 



with 



B\0) = f dk 9_(k, j3) d\ r {k) d\ r (-k) + 9 + (k, j3) d Ir (k) d Ir (-k) 



(14) 



The underlined operators are the operators corresponding to the doubling of the Hilbert 
space that arises in thermofield dynamics method. We shall determine the condensate 
function f(k), and the functions 9_(k,(3) and 9 + (k,[3) of the thermal vacuum through min- 
imisation of the thermodynamic potential. To evaluate the expectation value of the energy 
density with respect to the thermal vacuum, we shall use the following formula. 



(^(x)My))fs = -^x-J (A_(£,/?))^ e -*<*-*>d£, 



(15) 



where 



A_(fc,/J) 



fcos 9+ + sin 



1 °cos( X (k)-2f(k)) 



+ a ■ k sm( X (k) - 2f(k)) (cos 2 9 + - sin 2 9 



(16) 



We now proceed to calculate the energy density, 



e = {H)p = e N + e a + e u 



(17) 



with 



CiV 



7 



(2*) J 



dk 



e(Jb)cos2/(Jb)-^(Mcos2/(A-) + \k\ sin 2. f(k ) 



(cos 2 0+ - sin 2 9_ 



(18a) 



^ = 2 m l a li 



18b) 



and 



^w = 9u)^o 7(27r) / dk(cos 9 + + sin 9 



\ m l^l- 



:i8c) 



The thermodynamic potential is then given as 



tt = e- -S - fip B , 



(19) 



with the entropy density 

sin 2 6L ln(sin 2 6L) + cos 2 6L ln(cos 2 6> 



S a + S^ (20) 



5 = - 7 (27t)~ 3 fdk 

+ sin 2 6> + ln(sin 2 0+) + cos 2 6> + ln(cos 2 H 

and the baryon density 

p B = 7 (2tt)- 3 f dk(cos 2 0+ + sin 2 0_). (21) 

In the above, 7 is the spin isospin degeneracy factor and is equal to 4 for nuclear matter. 
Further, S a and S u are the contributions to the entropy density from a- and u;-mesons, 
respectively. It may be noted here that these are independent of the functions f(k), 9±(k, 13) 
associated with the nucleons and hence are not relevant for the nuclear matter properties at 
zero temperature. Extremising the thermodynamic potential Q with respect to the conden- 
sate function f(k) and the functions 9^ yields 

— t 

tan2/(fl= °°lf) (22) 

e{k) 2 + Mg a a 

and 

sin2 e T = 7777 /,v tt 7 (23) 

with e*(k) = (k 2 + M* 2 ) 1 / 2 and /i* = /j, — g w UQ as the effective energy density and effective 
chemical potential, where the effective nucleon mass is M* = M + g a a . 
Then the expression for the energy density becomes 

e = e N + e a + e w , (24) 

with 

/(2tt)- 3 j dk(k 2 + M* 2 ) 1 / 2 (sin 2 9_ - cos 2 0+), (25a) 

^ = 2 m ^o. ( 25b ) 



e^v = 71 



and 



r _ I 

^ = g u oj G 7(2vr) -3 / dfc(sin 2 9- + cos 2 9+) - -m 2 ^. (25c) 

We now proceed to study the properties of nuclear matter at zero temperature. In that 
limit the distribution functions for the baryons and antibaryons are given as 

sin 2 ^ = Q(i2* -e*(k)); sin 2 fl + = 0. (26) 

The energy density after subtracting out the pure vacuum contribution then becomes 

e„ = e(0_,/)-e(0_ = O,/ = O) 

= e MF T + Ae (27) 



with 



e M FT = 7(2tt)- 3 J d^h 2 + M* 2 ) 1 ' 2 + -mlol + g u u oPB - -ml^l (28) 



|fc|<fe F 



and 



Ae = - 7 (2vr)- 3 J dk\(k 2 + M* 2 ) 1 / 2 - (k 2 + M 2 ) 1 ' 2 



g a o- M 



(29) 



(P + M 2 ) 1/2 - 

The above expression for the energy density is divergent. It is renormalised by adding 
the counter terms 

eat = E °n°l (30) 

71=1 

The addition of the counter term linear in o"o amounts to normal ordering of the scalar 
density in the perturbative vacuum and cancels exactly with the last term in equation ( |29| ) 
0. The first two terms of the same equation corresponds to the shift in the Dirac sea arising 
from the change in the nucleon mass at finite density when a acquires a vacuum expectation 
value, and consequent divergences cancel with the counter terms of ([31]) with higher powers 
in o"o 0. Then we have the expression for the finite renormalised energy density 

tren = ^MFT + Ae ren , (31) 

where 



Ae r 



7 



16vr 2 



,M* 



*\2 



M* 4 ln( — ) + ikP(M - M*) - -M\M - M*) 
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+ —M(M - M*f M(M - M*) 4 



(32) 



For a given baryon density as given by 



p B = 1 (2n)- 3 / dkQ(k F -k) 



(33) 

the thermodynamic potential given by equation (plj) is now finite and is a function of <To 
and uq. This when minimised with respect to <jq gives the self-consistency condition for the 
effective nucleon mass, 

M* 



M* = M 



al 7 

ml (2tt) 



dk 



e(k)* 



Q(k F -k) + AM* 



(34) 



where 

AM* 



9a 7 

m 2 (2tt) 



,M* 
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M* 3 ln(— ) + M 2 (M- AT) - -M 2 (M - M*f + —M(M - M* v '' 



(35) 



We note that the self-consistency condition for the effective nucleon mass as well as 
the energy density as obtained here through an explicit construct of a state with baryon- 
antibaryon condensates are identical to those obtained through summing tadpole diagrams 
for the baryon propagator in the relativistic Hartree approximation 0. 

III. ANSATZ STATE WITH BARYON ANTIBARYON AND SIGMA 

CONDENSATES 



We next consider the quantum corrections due to the scalar mesons as arising from 
a vacuum realignment with sigma condensates. This means that the a-field is no longer 
classical, but is now treated as a quantum field. As will be seen later, a quartic term in 
the sigma field would favour such condensates. Self-interactions of scalar fields with cubic 



and quartic terms have been considered earlier [18| in the no- sea approximation as well 
as including the quantum corrections arising from the sigma fields p,19,20fl. They may be 
regarded as mediating three- and four-body interactions between the nucleons. The best 



9 



fits to incompressibility in nuclear matter, single-particle spectra and properties of deformed 
nuclei are achieved with a negative value for the quartic coupling in the sigma field. However, 
with such a negative coupling the energy spectrum of the theory becomes unbounded from 
below |21[ for large o and hence it is impossible to study excited spectra or to include vacuum 
polarisation effects. 

Including a quartic scalar self-interaction, eq. (|5b|) is modified to 

H a = ^ad»a + ^m 2 a a 2 + \a\ (36) 

with ma and A being the bare mass and coupling constant respectively. The a field satisfies 
the quantum algebra 

[<j{x), <y(y)] = iS(x - y). (37) 

We may expand the field operators in terms of creation and annihilation operators at time 
t = as 

a(x : 0) = — |Lj J -JL= (a(k) + a\-k)) e ** (38a) 



( 2 -) 3/2 J sj^$) 



u(k) 



a(x,0) = -^j- 2 Jdk^ ^Y (-a(k) + a\-k)) e^. (38b) 



In the above, io(k) is an arbitrary function which for free fields is given by uo(k) = \Jk 2 + va 2 a 
and the perturbative vacuum is defined corresponding to this basis through a \ vac) = 0. 
The expansions ( |3"8"D and the quantum algebra ([37]) yield the commutation relation for the 
operators a as 

\a{k),a\k')} = 5(k-k'). (39) 

As seen in the previous section a realignment of the ground state from | vac) to | vac') 
with nucleon condensates amounts to including quantum effects. We shall adopt a similar 
procedure now to calculate the quantum corrections arising from the cr-field. We thus modify 
the ansatz for the trial ground state as given by flTOl) to include a condensates as |13| 
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\Sl) = U U F \vac), 



(40) 



with 



U 9 = U n Ui 



where Ui = exp(BJ — B t ), (i = I, II). Explicitly the B { are given as 



(41) 



B\ 



dk 



\ 



^-UkW(k), 



(42a) 



and 



B H ] = - J dkg(k)a'\k)a'\-k). 



(42b) 



u(k) 



In the above, a'(k) = Uia{k)Uj = a(k) — \i2^fa(k) corresponds to a shifted field operator 
associated with the coherent state [13|] and satisfies the same quantum algebra as given in 
eq. (0). Thus in this construct for the ground state we have two functions f a (k) and 
g(k) which will be determined through minimisation of energy density. Further, since | Q) 
contains an arbitrary number of a'' quanta, a' \ O) ^ 0. However, we can define the basis 
b(k), b*(k) corresponding to | Q) through the Bogoliubov transformation as 



' b(k) ^ 



b\-k) 



U, 



1 a' $) ^ 



n 



\ 



-k) 



UJi 1 



( 
\ 



cosh(yf — smhg 
- sinhg coshg 



\ ( a'(k) ^ 



V 



■k) 



(43) 



It is easy to check that b(k) \ Q) = 0. Further, to preserve translational invariance f a (k) 
has to be proportional to 8(k) and we take f a (k) = ao(2ir) 3 / 2 5(k). o"o will correspond to a 
classical field of the conventional approach ||13| . We next calculate the expectation value of 
the Hamiltonian density for the cr-meson given by equation(^B|). Using the transformations 
(fi3D it is easy to evaluate that 



(Q \ a \ Q) = o"o, 



(44a) 



but, 



(SI | a 2 | fi) = a 2 + I, 



(44b) 
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where 

1 f dk 
/= (2^/2-^) (c ° sh29 + sinh29) - (44C) 

Using equations (|36|) and (44) the energy density of H a with respect to the trial state 
becomes [|13 



11 r dk 



^<n|«.|n> = 27£i*/ 



2 (2vr) 3 J 2u{k) 



k 2 (smh.2g + cosh2g) + u 2 (k)(cosh.2g — sinh2g) 



+ -mil + QXa 2 I + 3A/ 2 + -m 2 a a 2 + Aa*. (45) 

Zi Zi 

Extremising the above energy density with respect to the function g(k) yields 

, ^ /, n 6AJ + 6Ao- 2 /^x 

tanh.2g(k) = , . 46 

uo{kf + 6AJ + 6A<x 2 
It is clear from the above equation that in the absence of a quartic coupling no such con- 
densates are favoured since the condensate function vanishes for A = 0. Now substituting 
this value of g(k) in the expression for the a- meson energy density yields 

e a = \mW + W + \j^j dk{k 2 + Mlfl 2 - 3XI 2 (47) 



where 



M 2 a =ml + 12X1 + l2Xa 2 (48) 



with 



1 fdk 



(49) 



(2tt) 3 7 2 (P + M2)V2 
obtained from equation (|44cD after substituting for the condensate function g(k) as in equa- 



tion (|46|) . The expression for the "effective potential" e CT contains divergent integrals. Since 
our approximation is nonperturbatively self-consistent, the field-dependent effective mass 
M a is also not well defined because of the infinities in the integral / given by equation (|49"|). 
Therefore we first obtain a well-defined finite expression for M a by renormalisation. We use 



the renormalisation prescription of ref. |£2[ and thus obtain the renormalised mass m# and 
coupling Xr through 

12 



^ = ^ + 12/! (A), (50a) 



^- = \ + 12I 2 {A,fj), (50b) 

ar A 



where I± and I2 are the integrals 



IfcKA 



|fc|<A 

with jjl as the renormalisation scale and A as an ultraviolet momentum cut-off. It may be 
noted here that with the use of the above renormalisation prescription the effective sigma 
mass M a and the energy density ultimately become independent of A and stay finite in the 
limit A — > 00. Using equations ( |50a| ) and (|50b| ) in equation (f48[) , we have the gap equation 
for Ml in terms of the renormalised parameters as 

Ml = m 2 R + \2\ R al + l2X R I f (M a ), (52) 

where 

Then using the above equations we simplify equation (|47 ) to obtain the energy density for 
the a in terms of <To as 

_3A^ + I g-)%^( ln (f)-I)-3 M - 2 K. (54) 

The above expression is given in terms of the renormalized a mass tur and the renormalized 
coupling A/? except for the last term which is still in terms of the bare coupling constant A 



and did not get renormalised because of the structure of the gap equation p3| . However, 
from the renormalisation condition (|50b| ) it is easy to see that when X R is kept fixed, as the 
ultraviolet cut-off A in eq. ( |51b|) goes to infinity, the bare coupling A —> 0_. Therefore the 
last term in eq. fl54|) will be neglected in the numerical calculations. 
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After subtracting the vacuum contribution, we get 
Ae CT = e CT - e CT (a = 0) 



1 2 2j o, *± M U\ ( M °\ l \ *\ t2 
= z m R a o + 3A fi a + — (In (—J - -J - 3A*!, 

where M CT)0 and Ijq are the expressions as given by eqs. (|52|) and (]53|) with cr = 0. 

In the limit of the coupling, Xr = 0, one can see that eq. (|55|) reduces to eq. fl25b| ) as it 



should. Also, we note that the sign of Xr must be chosen to be positive, because otherwise 
the energy density would become unbounded from below with vacuum fluctuations [B,20,21 . 



The expectation value for the energy density after subtracting out the vacuum contribu- 
tion as given by eq. (P7|), now with sigma condensates is modified to 

eo = e { mite + Ae, (56) 

where 

6 /mite = 7 (2vr)- 3 j dk(k 2 + M*Y 2 + 9^oPB-lmlu; 2 +Ae a , (57) 

\k\<k F 

with Ae CT given through eq. ([55]) and Ae is the divergent part of the energy density given by 
equation (p9|). We renormalise by adding the same counter terms as given by (^) so that 
as earlier the renormalised mass and the renormalised quartic coupling remain unchanged 
0. This yields the expression for the energy density 

with Ae ren given by eq. (|32j). As earlier the energy density is to be minimised with respect 
to o" to obtain the optimised value for cr , thus determining the effective mass M* in a 
self-consistent manner. 

The energy density from the a field as given by eq. (|55|) is still in terms of the renormali- 
sation scale /i which is arbitrary. We choose this to be equal to the renormalised sigma mass 
tur in doing the numerical calculations. This is because changing /i would mean chang- 
ing the quartic coupling Xr, and Xr here enters as a parameter to be chosen to give the 
incompressibility in the correct range. 
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For a given baryon density, ps, the binding energy for nuclear matter is 

E B = e ren /p B - M, (59) 

The parameters g a , g^ and Xr are fitted so as to describe the ground-state properties of 
nuclear matter correctly. We discuss the results in the next section. 

IV. RESULTS AND DISCUSSIONS 

We now proceed with the numerical calculations to study the nuclear matter properties 
at zero temperature. We take the nucleon and o;-meson masses to be their experimental 
values as 939 MeV and 783 MeV. We first calculate the binding energy per nucleon as 
given in equation ( |59| ) and fit the scalar and vector couplings g a and g u to get the correct 
saturation properties of nuclear matter. This involves first minimising the energy density in 
eq. (|5^) with respect to a to get the optimised scalar field ground state expectation value 
u min . This procedure also naturally includes obtaining the in- medium a- meson mass M a 
through solving the gap equation (p2|) in a self-consistent manner. Obtaining the optimized 
a min amounts to getting the effective nucleon mass M* = M + g a a min . We fix the meson 
couplings from the saturation properties of the nuclear matter for given renormalised sigma 
mass and coupling m_R and \r. Taking rriR = 520 MeV, the values of g a and g u are 7.34 
and 8.21 for X R = 1.8, and are 6.67 and 7.08 for A^ = 5 respectively. Using these values, we 
calculate the binding energy for nuclear matter as a function of the Fermi momentum and 
plot it in fig. 1. In the same figure we also plot the results for the relativistic Hartree and 
for the no-sea approximation. Clearly, including baryon and cr-meson quantum corrections 
leads to a softer equation of state and the softening increases for higher values of \r. The 



incompressibility of the nuclear matter is given as 121] 



evaluated at the saturation Fermi momentum. The value of K is found to be 401 MeV for 
\r = 1.8 and 329 MeV for \r = 5. These are smaller than the mean-field result of 545 MeV 
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], as well as that of relativistic Hartree of 450 MeV and are similar to those obtained 



in ref. |2(J containing cubic and quartic self-interaction of the a-meson. 

In fig. 2, we plot the effective nucleon mass M* = M + g a <J m in as a function of Fermi mo- 
mentum with a m in obtained from the minimisation of the energy density in a self-consistent 
manner. At the saturation density of k F = 1.42 fin -1 , we get M* = 0.752M and 0.815M 
for \r = 1.8 and \r = 5, respectively. These values may be compared with the results of 
M* = 0.56M in the no-sea approximation and of 0.72M in the relativistic Hartree. 

In fig. 3, we plot the vector and the scalar potentials as functions of k F for sigma self 
coupling Ar = 1.8 and 5. At saturation density the scalar and vector contributions are 
U s = g a a min = -232.7 MeV and U v = g u u = 163.4 MeV for X R = 1.8 and are -173.14 
MeV and 107.74 MeV for \r = 5 respectively. These give rise to the nucleon potential 
(Us + Uy) of —69.3 MeV and —65.4 MeV and an antinucleon potential (Us — Uv) of —396.1 
MeV and —280.9 MeV for Ar = 1.8 and 5. Clearly the inclusion of the quantum corrections 
reduces the antinucleon potential as compared to both the relativistic Hartree (—450 MeV) 
and the no-sea results (-746 MeV) §. 

In fig. 4, we plot the in-medium cr-meson mass M a of eq. ( [52] ) as a function of baryon 
density for X R = 1.8 and 5. M a increases with density as \r is positive and the magnitude 
of <J m i n increases with density too. However, the change in M a is rather small. 

In fig. 5, we plot the incompressibility if as a function of the quartic coupling X R for 
different values of ttir, the renormalised sigma mass in vacuum. The value of K decreases 



with increase in \ R similar to the results obtained in ref. [ 20f| . In fig. 6, we plot the effective 



nucleon mass versus the sigma self coupling for various values of tur. The value of M* 
increases with Xr, which is a reflection of the diminishing nucleon-sigma coupling strength 
for larger values of the quartic self-interaction. 

To summarise we have used a nonperturbative approach to include quantum effects in 
nuclear matter using the framework of QHD. Instead of going through a loop expansion and 
summing over an infinite series of Feynman diagrams we have included the quantum correc- 
tions through a realignment of the ground state with baryon as well as meson condensates. 
It is interesting to note that inclusion of baryon-antibaryon condensates with the particular 
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ansatz determined through minimisation of the thermodynamic potential yields the same 
results as obtained in the relativistic Hartree approximation. This results in a softer equa- 
tion of state as compared to the no-sea approximation. Inclusion of scalar meson quantum 
corrections in a self-consistent manner leads to a further softening of the equation of state. 
The value for the incompressibility of nuclear matter is within the range of 200-400 MeV 



25fl . It is known that most of the parameter sets which explain the ground state properties 
of nuclear matter and finite nuclei quite well are with a negative quartic coupling. But the 
energy spectrum in such a case is unbounded from below [21 1 for large a thus making it 



impossible to include vacuum polarisation effects. We have included the quantum effects 
with a quartic self interaction through sigma condensates taking the coupling to be positive. 
We have also calculated the effective mass of the sigma field as modified by the quantum 
corrections from baryon and sigma fields. The effective sigma mass is seen to increase with 
density. 

We have also looked at the behaviour of the incompressibility as a function of the coupling 
Xfi for various values of sigma mass, which is seen to decrease with the coupling. Finally, we 
have looked at the effect of the sigma quartic coupling on the effective nucleon mass which 
grows with the coupling. Generally, higher values of the quartic term in the potential of 
the cr-meson tend to reduce the large meson fields and thus the strong relativistic effects 
in the nucleon sector. Clearly, the approximation here lies in the specific ansatz for the 
ground-state structure. However, a systematic inclusion of more general condensates than 
the pairing one as used here might be an improvement over the present one. The method 
can also be generalised to finite temperature as well as to finite nuclei, e.g., using the local 
density approximation. Work in this direction is in progress. 
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Figures Captions 

Fig. 1. The binding energy for nuclear matter as a function of Fermi momentum kp cor- 
responding to the no-sea and the relativistic Hartree approximations and the approach 
including quantum corrections from baryon and cr-meson given by eq. (59). It is seen that 
the equation of state is softer with such quantum corrections. 

Fig. 2. The effective nucleon mass for nuclear matter as a function of the Fermi momentum, 

Fig. 3. The scalar potential Us (negative values) and vector potential Uy (positive values) 
for nuclear matter as functions of Fermi momentum kp. 

Fig. 4. The in-medium cr-meson mass M a of eq. (|52|) as a function of density, which is seen 
to increase with density. However, the change is seen to be rather small. 

Fig. 5. The incompressibility K versus the quartic coupling \r for various values of m#, 
which is seen to decrease with A#. The values of K are higher for larger values of the sigma 
mass. 

Fig. 6. The effective nucleon mass as a function of \r for different values of tur. It is seen 
to decrease with increase in the coupling. 
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Figure 2 
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Figure 3 
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Figure 4 
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Figure 5 
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